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ABSTRACT 

Models for the formation and evolution of galaxies readily predict physical properties such as the star for- 
mation rates, metal enrichment histories, and, increasingly, gas and dust content of synthetic galaxies. Such 
predictions are frequently compared to the spectral energy distributions of observed galaxies via the stellar pop- 
ulation synthesis (SPS) technique. Substantial uncertainties in SPS exist, and yet their relevance to the task of 
comparing galaxy evolution models to observations has received little attention. In the present work we begin 
to address this issue by investigating the importance of uncertainties in stellar evolution, the initial stellar mass 
function (IMF), and dust and interstellar medium (ISM) properties on the translation from models to observa- 
tions. We demonstrate that these uncertainties translate into substantial uncertainties in the ultraviolet, optical, 
and near-infrared colors of synthetic galaxies. Aspects that carry significant uncertainties include the loga- 
rithmic slope of the IMF above IMq, dust attenuation law, molecular cloud disruption timescale, dumpiness 
of the ISM, fraction of unobscured starlight, and treatment of advanced stages of stellar evolution including 
blue stragglers, the horizontal branch, and the thermally-pulsating asymptotic giant branch. The interpretation 
of the resulting uncertainties in the derived colors is highly non-trivial because many of the uncertainties are 
likely systematic, and possibly correlated with the physical properties of galaxies. We therefore urge caution 
when comparing models to observations. 

Subject headings: galaxies: evolution — galaxies: stellar content 



1. INTRODUCTION 

Models for the formation and evolution of galaxies have be- 
come dramatically more sophisticated over the past decade. 
These models span a variety of techniques, including the 
semi-analytic approach, which couples dark matter halo 
merger trees to analytic recipes for the evolution of th e 
baryons (e.g. IWhite & Frenkl 1 1 99 U iKauffmann eTall 1 1 9991) . 
and cosmologically-embedded hydrodynamic simulations, 
which numerically follow the evolution of dark matter and 
baryons self-consistently, with numerical recipes encompass- 
ing unresolved physical processes such as star formation, 
black hole growth, stellar and black hole feedback, etc (e.g. 
iCen & Ostrikerlll992HKatz et al.lfl99l iKeres et al.ll2005h . 

The primary predictions of these models include quan- 
tities such as the total mass in stars, star formation and 
metal enrichment histories, and, in some cases, informa- 
tion on cold gas properties. Unfortunately, these quanti- 
ties are not directly available to extragalactic observers. In- 
stead, these physical quantities are often transformed into 
observable spectral energy distribution s via stellar popula- 
tion syn t hesis ( SPS) techniques (e .g. Tinslev & Gunn! [1976t 
Bruzu al 119831: iRenzini & Buzzon f 1 19861 Worthevrll994 
Fioc & Rocca-Volmerange 1997; Maraston 1998; Vazdekis 
1999). For many applications, SPS thus provides the bridge 



Modeling of the attenuation of starlight by dust is a partic- 
ularly challenging aspect in the SPS technique. Such mod- 
els fall into roughly two categories. On the one hand there 
are approac hes that use physical models for the properties o f 
grains (e.g. IWeingartner & Draindl2001t IZubko et afl 120041) . 
coupled to assumptions ab out the geometry of the dust with 
respe c t to the stars (e.g. ISilva et all 1 1998t Ipevriendt et al.1 
2001; T uffs et al.l 12001; |Jonssodl2006f 
These models rely on radiative trans- 



1999; Gordon et al. 
Galliano_elalJ[2008) 



between models and observations. 

Until now, little attention has been paid to the uncertain- 
ties in the SPS technique in the context of translating galaxy 
formation model predictions into observables. And yet uncer- 
tainties in advanced stages of stellar evolution, stellar spec- 
tral libraries, the stellar initial mass function, and the adopted 
model for the obscuring effects of interstellar dust can each 
be significant components of the total error budget on the pre- 
dictions. 



fer calculations and are therefore computationally expensive. 
Their computational cost has given rise to the second class of 
models, which are physically-motivated, but are in essence 
phenomenological. Models of the second class typically 
adopt an attenuation curve with a normalization, and, in prin- 
ciple, shape, that d epends on the age of the system (e.g. 
ICharlot & F all 2000). This effective attenuation curve is then 
applied uniformly to all stars of a given age. 

Common to both approaches is the fact that the required in- 
puts are in many cases seriously underconstrained. An accu- 
rate dust model requires knowledge of the grain composition 
and size distrib ution, and the wavelength-dependence of the 
grain albedo (Draine 2003), yet these quantities are extremely 
diffic ult to constrain in the Galaxy (e.g. Weingartner & Draine 
120011) . let alone in external systems. It is expected that grain 
properties will be functions of metallicity and the local inten- 
sity of ultraviolet radiation — in addition to other variables — 
both of which will vary from galaxy to galaxy. The geome- 
try of the dust, including large-scale inhomogeneities and the 
proximity of the dust to the stars, is not known at all outside of 
the Local Group, and not known with the requisite precision 
in our own Galaxy. Yet the geometry of dust has a substan- 
tial effect on the net attenuation of starlight. The lifetime of 
molecular clouds is another important parameter because it 
controls the length of time that young stars are heavily dust 
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obscured, and yet this timescale is not known to better than 
an order of magnitude, and will likely depend on quantities 
such as metallicity and local st ar formation rate (Blitz & Shul 
1980 1- iMcKee & Ostrikei]|2007|). 

iFontanot et al.l (120091) have recently compared a physical 
dust model that employs radiative transfer to simple phe- 
nomenological prescriptions. These authors find substantial 
differences between the analytic calculations and the more 
complex physical dust model. While comparisons of this type 
are essential for informing analytic prescriptions, our belief is 
that the physical dust models themselves are sufficiently un- 
certain to warrant a flexible approach to the obscuring effects 
of dust. 

Uncertainties in stellar evolution calculations can also im- 
pact t he translation from models to obser yables (e.g. iTinslev I 
M980) ICharibt et all 1 19961: ICharlotl ll996l iLee et al.ll2007t lYil 
l2003l:IConroy et al.ll2009l) . Recently lTonini et al.l (120091) have 
investigated the impact of the thermally-pulsating asymptotic 
giant branch (TP-AGB) phase on the observational predic- 
tions of a semi-analytic galaxy evolution model. They com- 
pare a popular SPS model that does not incl ude the TP-AGB 
phase of stellar evolution to the model of Mara stonl (120051) 
where this important but uncertain phase is handled with care. 
These authors find that the near-IR colors of model galaxies 
can differ by as much as two magnitudes between the two SPS 
models, and that the differences increase with redshift. These 
results highlight the importance of carefully accounting for 
uncertain aspects of SPS when translating galaxy formation 
model predictions into observables. 

In the present work we extend the spirit of Tonini et al. by 
considering a broader array of uncertainties in stellar evolu- 
tion, uncertainties in the stellar initial mass function (IMF), 
and the substantial uncertainties associated with the treat- 
ment of interstellar dust. We focus on two synthetic galaxies 
drawn from a recent semi-analytic model of galaxy evolution. 
These galaxies were chosen to represent a typical bright star- 
forming galaxy and a typical bright passively evolving galaxy, 
both at z = 0. For the purposes of this work, the main differ- 
ence between these two types of galaxies is the presence or 
absence of young, hot stars and the average stellar age. 

This paper is the second in a series focusing on uncer- 
tainties in stellar population synthesis. In this series, par- 
ticular attention is focused on how those uncertainties effect 
o ur knowledge of ga laxy formation and evolution. In Paper 
I dConrov et alj 120091) we examined the impact of SPS un- 
certainties on deriving physical properties of observed galax- 
ies, including star formation rates, ages, and stellar masses. 
This paper considers the inverse problem: the uncertainties in 
translating model galaxies into observables. 

This paper continues with $2] where we discuss the salient 
ingredients of SPS, including a detailed description of our 
dust treatment, ^briefly describes the synthetic galaxies and 
the semi-analytic model from which they are derived. Sj4]con- 
tains our main results, where we investigate the uncertainties 
in the predicted UV, optical, and near-IR colors of synthetic 
galaxies due to uncertainties in SPS. Following this sobering 
assessment, in S|5]we comment on several implications. A 
summ ary is provided in $6] All magnitudes are in the AB sys- 
tem (lOke & Gunnll983l) . Where necessary a Hubble constant 
of Ho = 100/; km s -1 Mpc" 1 is assumed. 

2. STELLAR POPULATION SYNTHESIS 

Our SPS treatment closely follows that of IConrov et al.l 
(2009), to which the reader is referred for details. In 



brief, the SPS code u ses the latest stellar evolution tracks 
from the Padova group (Marigo & Girardi 2007; Marig cTet al.l 
2008), which follow stellar evolution from the main se- 
quence through the thermally-pulsating asymptotic giant 
branch (TP-AGB) phase. Evolutionary calculations exist for 
metallicities in the range 10" 4 < Z < 0.030, for ages 10 6 ' 6 < 
t < 10 10 2 yrs, and for initial masses 0.15 < M < 67 M . The 
stellar spectral libraries are primarily thos e of the empirically- 
calibrated theoretical BaS eL3.1 library dLeieune et al.lll997l 
Il998tlw"estera et aT] |2002), supp lemented with empirical TP- 
AGB spectra from the library of Lancon & Mouhcine (2002). 
The initial stellar mass function (IMF) of KroupaT ikOQll) is 
adopted. Our SPS code is open-source and publicly avail- 
able 1 . 

We seek to assess the importance of uncertainties in SPS 
in translating the predictions of galaxy evolution models into 
observables. In our approach we will focus on two distinct 
classes of uncertainties. The first class are uncertainties in 
the single stellar populations (SSP) resulting from isochrone 
synthesis, including both the stellar evolution tracks and the 
IMF. The second class concerns the treatment of reddening by 
interstellar dust. 

2.1. SSP uncertainties 
2.1.1. Stellar evolution ( isochrone ) uncertainties 

Uncertainties in the isochrones are quantified as in 
IConrov et al1d2009l) . The position of TP-AGB stars in the HR 
diagram is substantially uncertain, both observationally and 
theoretically. We thus introduce two variables that amount to 
shifts in log(r e ff) and log(Lboi) with respect to the default stel- 
la r evolution tracks: A j- and Al, respectively. As discussed 
in lConrov et al.1 (120091) . these variables effectively encompass 
uncertainties not only in the stellar evolution tracks but also 
uncertainties in the associated spectral energy distributions 
and circumstellar dust enshrowding these stars. Since these 
stars have relatively low effective temperatures, they con- 
tribute substantially to integrated spectra only at A > 7000A. 

Blue straggler (BS) stars and blue horizontal branch (BHB) 
stars are almost universally neglected when translating galaxy 
evolution models into observables, although their importance 
has been demonstrated w h en modeling obs e rved spectra (see 
e.g. iJimenez et all 120041; iMarastonl l200l iLi & Hani 120081. 
2009; Con rov et al.l 120091) . We introduce two additional pa- 
rameters to quantify these uncertain stellar phases: the spe- 
cific frequency of BSs, 5bs, defined as the number of BSs per 
unit HB star, and the fraction of HB stars that are blue, /bhb- 
In this context a blue HB star is any HB star extending blue- 
ward of the red clump identified in standard stellar evolution 
tracks. In our treatment this extended blue component is pop- 
ulated uniformly in log(r e ff) from the red clump to T e s = 10 4 
K. 

We do not allow for any of our parameters to be age- 
dependent, but note that one may interpret a given parame- 
ter value as being representative of some effective population 
age. 

2.1.2. IMF uncertainties 

We will also explore the importance of the logarithmic 
slope of the IMF on the derived colors of galaxies. It is widely 
understood that the IMF has a significant impact on the mass- 
to-light ratio of galaxies, but it's impact on the SEDs of galax- 

1 www.astro.princeton.edu/'^cconroy/SPS/ 
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ies is less appreciate d (or rather, h as become less appreciated 
with time, see e.g., Tinslevl [l980l) . It is important to consider 
the influence of the IMF on the translation between models 
and observations because of the difficulty in measuring the 
IMF in the sola r neighborhood , let alone in external galaxies. 

For example, iRroupal (1200 ll) compiled a number of obser- 
vations of the logarithmic slope of the IMF in the solar neigh- 
borhood and found a = 2.3 ± 0.7 for m > IM©, at 99% con- 
fidence, where the IMF is parameterized as £(m) oc m~ a . It is 
important to realize that th is uncertainty is only statistical. As 
discussed by Kroupa (2001), corrections due to unseen binary 
companions will depend on a, and will result in an increase 
in a of the order of 0.3, depending on the details of the binary 
companions. Further complications arise near IM© because 
the observational techniques for determining the IMF at this 
mass regime rely on star counts in the field. At ~ IM© the 
stellar lifetime is of order the age of the Galactic disk, and 
so evolutionary corrections, which are substantial, rely on a 
detailed knowledge of the star formation history of the disk. 
In addition, because the age of IM© stars is of order a Hub- 
ble time, the IMF near this mass matters for understanding 
galaxy SEDs. The uncertainty on a at ~ 1M© is so lar ge that 
data p oints near 1M© are often excluded from fits (Krourja 
12001 . 

The physical conditions of the Universe relevant to star for- 
mation were ve ry different at h igher redshift, including higher 
ISM p ressures dLiu et al.l l2008). lower metallicities dErb et al.l 
2006), and higher CMB temperatures. The latter will start 
to effect the f ragmentation and collapse of gas around z ~ 5 
(lLarsonll2005l) . Some authors have suggested that these con- 
ditions will result in a top-heavy (or bottom-light) IMF at 
earlier times, owing essentially to a larger Jeans mass (e.g. 
Larson 1998, 2005). Recent observational evidence, includ- 
ing th e abundance patterns of metal-poor s tars in the Milky 
Way dLucatello et al.l2005HTumlinsonl20 07). evolution in the 
fundamental plane (van Dokkum 2008), and a comparison be- 
tween the cosmic mass density evo lution and s tar formation 
rate density evolution beyond z ~ 2 (Dave 2008), is beginning 
to support this picture, however tentatively ( for an alternative, 
less e xotic explanation of the last point, see Reddv & Steidel 
l2009h . 

Recent observational results in the local Universe has 
also pointed toward an IMF that varies with environment. 
Meur eFet al.l (|2009) recently demonstrated that the ratio be- 
tween Ha flux and near-UV luminosity varies systematically 
with galaxy surface brightness. This ratio should be constant 
for a universal IMF. The data can be explained by an IMF 
slope at high masses that varies over the range 1 .5 < a < 4.0, 
or by a systematically varying upper mass cut-off over the 
range w 20 < M up < 100M©. 

We follow Kroupa (2001) in the parameterization and best- 
fit logarithmic slopes of the IMF. In the following sections we 
will explore the effects of varying the IMF only for masses 
> IM©. Variation is limited to > IM© because the main- 
sequence turn-off mass is IM© for a 12 Gyr SSP, and so stars 
of lower mass never contribute substantially to the integrated 
colors (unless pathologi cal IM Fs are considered). In addition, 
as will be discussed in £14.1.21 it is primarily galaxies with a 
range of stellar populations that are sensitive to IMF varia- 
tions. In such galaxies it is the relative frequencies of stars 
with mass > IM© that contribute to the integrated light and 
therefore only IMF variations at > IM© will have impact on 
the colors. 



2.2. Dust treatment 

We now turn to the treatment of dust in the SPS technique. 
Since we are interested in restframe wavelengths bluer than 
w 2.4/im (i.e. the Zf-band), we do not model the reprocessing 
of radiation by dust; we only attempt to model its obscuring 
effects. 

The light emerging from a galaxy can be modeled as 



F x (t)= [ nt-t^Sxit'.Zit-t'^e-^At', 

JQ 



(1) 



where ^ is the SFR, S\ is the SSP spectrum, F\ is the result- 
ing composite spectrum, t\ is the optical depth, and the inte- 
gration variable is the age of stellar populations. Note that this 
formulation allows for metallicity evolution with time, Z(f )■ 

The obscuring effects of dust are encompassed in the opti- 
cal depth T\(t). A toy model for the effective optical depth 
was presented in lCharlot & Falll (120 00) and has since become 
quite popular 2 . The physical motivation for their model is 
the following. Stars are born in molecular clouds and thus 
the light from young stars is heavily attenuated by dust in the 
cloud. At later times the stars will no longer experience atten- 
uation due to their birth cloud, either because they will have 
wandered out of the cloud, or because the cloud will have 
evaporated. Stars are thus subject to attenuation that varies 
with time — young stars are obscured by the dense molecu- 
lar clouds in which they form in addition to the diffuse ISM, 
while the light from older stars is attenuated only by the dif- 
fuse ISM. Motivated by this physical picture, the optical depth 
is parameterized as: 



. n(X/ 5500k)- 5 ' 
' A S r 2 (A/5500A)-^ 



t < t a 
t >f es 



(2) 



where f esc is the timescale over which young stars reside in 
their natal clouds. While there is ample evidence t hat the 
extinc tion curve slope depends on environment (e.g. Mathis 
1990), herein we fix J| = 82 = 5, which ha s become a com- 
mon a ssumption (e.g. ICharlot & Fallll2000T) 3 . ICharlot & Fal 
(2000) advocate 5 = 0.7, log(f esc /yrs) = 7.0, n = 1.0, and 
T2 = 0.3, although it is clear from their own figures that these 
parameters may plausibly vary over a wide range. We choose 
to leave these parameters free w ith a range motivated by a 
variety of data (see discussion in $2.3\ . 

In reality each of the parameters t\ , r%, S, and f esc will take 
on a range of values for a given galaxy. For example, t\ 
will depend on the location of the young star cluster with re- 
spect to the molecular cloud within which it's embedded. If 
the stars are born near the edge of the cloud, as is often ob- 
served within the Galaxy (llsraell l 1978). then t\ will depend 
strongly on the relative orientation of the observer with re- 
spect to the star cluster — cloud configuration. A more accu- 
rate dust model would therefore allow for a distribution in 
each of these quantities. Owing to the non-linear relation be- 
tween these parameters and the resulting galaxy spectrum, the 
use of distributions rather than a mean values will have impor- 
tant consequences. These issues will be considered in detail 



2 Th is is one of the dust models provided in the popular Bruzual & Chariot 
12003) SPS code, and therefore is a standard dust model for many investiga- 
tors. 

3 In future work we will investigate this assumption in detail. For the 
present work, we are interested in quantifying the uncertainties associated 
with the dust model, and so it matters little if we vary &\ and 82 separately or 
in conjunction. 
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in future work; in the present work we consider only mean 
values for t\ and <5. 

The optical depth associated with the diffuse ISM, T2, will 
be handled separately. Many treatments of dust attenuation 
adopt the reasonable assumption that the optical depth of 
dust is proportional to the column de nsity of metals (e.g. 
iGuiderdoni & Rocca-Volmerangd [l987). Of the parameters 
listed above, T2 is therefore the parameter most readily pre- 
dicted by galaxy evolution models. In our main results, we 
therefore choose not to marginalize over this important pa- 
rameter. It must be noted, however, that this parameter is not 
at all easily or reliably predictable in the majority of galaxy 
evolution models owing to the enormous dynamic range re- 
quired to adequately resolve the relevant physical processes 
leading to the production of metals. And, even if the quantity 
and spatial distribution of metals an be accurately predicted, 
the relationship between metals and dust is in reality quite 
complex owing to the variety of heating and cooling processes 
in the ISM relevant for dust formation and destruction. 

The attenuation law required in our formulation represents 
the net fraction of starlight photons removed by dust, as seen 
by a distant observer. It is different than true extinction, which 
quantifies the fraction of photons removed along the line of 
sight to a single star. The difference is that extinction is due 
to both true absorption and scattering, while attenuation is pri- 
marily true absorption and geometric effects. In other words, 
while photons from a single star may be scattered out of the 
line of sight to an observer, when viewing an entire galaxy 
those scattered photons will on average (depending on geome- 
try) be seen by a distant observer. Extinction can be measured 
in the Milky Way and Magellanic Clouds, but we require at- 
tenuation when applying dust corrections to external galaxies, 
and therefore require accurate know ledge of the net scattering 
effects of dust. See ICalzettil d2~00 1 ) for detailed discussion of 
these and related issues. 

It is well-known that the Milky Way and Magellanic Cloud 
extinction curves are approximately a one parameter family 
dCardelli et all 119891: lFitzpatrick|[T999l) . In the optical and 
near-IR these curves are well-characterized by a power-law. 
The net attenuation in local starburst galaxies is also approxi- 
mately a po wer-law from the far-UV to « 1 ym (ICalzetti et alJ 
1994, 2000). Motivated by the qualitative aspects of these ob- 
servational results, we have adopted a power-law attenuation 
law in our dust model. For the purposes of this paper it is not 
essential that our parameterization provide the most accurate 
reflection of the underlying attenuation curve; rather, it is im- 
portant only that the parameterization be flexible enough to 
encompass the uncertainties, or variations, in the attenuation 
law. Or, put differently, for our purposes we are primarily in- 
terested in an accurate parameterization of the variance, and 
not the mean, of the attenuation law. 

There is an additional feature in the Galactic and Magel- 
lanic Cloud extinction curves that deserves mention. The 
strong, broad absorption feature at 2 175 A is ubiquitous in 
the Galactic and LMC extinction curves, and has even been 
detected at z ~ 2 dNoll & Pierinill2005l;lEhasd6ttir et aT1l2009l; 
iNoll et alj2009l) . It is however conspicuo usly absent in the at- 
tenuation curves within local starbursts (ICalzetti et al.ll 1994b 
and in the extinc tion curves through most sightlines of the 
SMC (IPeH 119921) . It is thought th at this feature is due to 
polycyclic aromatic hydrocarbons dDraind [2003h . and that 
it may vary with lo cal environment and star formation rate 
dGordon et alj|2003h . 

This feature happens to lie within the GALEX NUV band 



for low-redshift galaxies. Its presence or absence can change 
the attenuation in the NUV band by a magnitude or more. We 
do not consider the uncertainties associated with this feature 
herein, but simply note that including it would likely double 
the resulting uncertainties in the NUV band. Future work will 
consider this issue in further detail. 

2.2.1. Large-scale dust distribution 

The prescription described above assumes that stellar pop- 
ulations are obscured by a uniform screen of dust. In real- 
ity, one expects the large-scale diffuse ISM to exhibit consid- 
erable variation. This effect is important because a clumpy 
ISM will result in a lower effective optical depth compared 
to a uniform screen, for the same total amount of gas and 
dust. Moreover, the effect of dumpiness is wavelength- 
dependent, implying that the dumpiness o f the ISM cannot 
be treated as a simple m odulation of T2 dNatta & P anagia 
Il984t IWitt & Gordon! 1 19961) . An accurate accounting of the 
large-scale inhomogeneities of the I SM is therefore a cru- 
cial element of any dust model (e.g. iBruzual A. et aTl l 1988 h 
I Witt et al.lll992l:IVarosi & Dwekll 1999 1: 1 Witt & Gordonl bood 
lGordonetal.ll200U Uonssonl 120061) . We now describe our 
treatment of a clumpy ISM. 

We define n as the density enhancement of the ISM along a 
given sightline, relative to a uniform screen. Then Equation!]] 
becomes: 



Fx(t) = 



P{n)e~ nTx[,l) f(f - t')S x (t',Z(t - 1 ')) dn df ', 



o Jo 



(3) 

where P(n) is the probability that an SSP is behind a column 
density of dust enhanced by a factor n relative to the uniform 
screen. The above relation holds as long as the optical depth 
is proportional to the column density of dust. 

We do not modify the treatment of obscuration around 
young stars where t < f esc . Therefore: 



Pin) i 



Sin- 
Pi 



1) t < f esc yr 

t > fesc yr, 



where S(x) is the Dirac delta function. 
There is a constraint that must be satisfied for Py. 



n P2 dn = 1 . 



(4) 



(5) 



This equation ensures that the clumping of dust does not alter 
the total amount of dust in the galaxy. 

The parameter T2 takes on a somewhat different meaning in 
the context of a clumpy ISM. Formally, it describes the atten- 
uation that would be experienced if the ISM were smoothed to 
a uniform density. Physically, T2 can be thought of as simply 
controlling the total dust mass in the galaxy — e.g., keeping 
T2 fixed while varying P2 implies that the total amount of dust 
is held constant, and only its spatial distribution is changed. 

There are very few direct observational constraints on 
the column density distribution o f gas and dust in the 
Galaxy , let alone in external systems. Berkhuiisen & Fletcher! 
(2008) measured the PDF of diffuse atomic gas in the 
solar neighborhood. They found that the distribution of 
densities is consistent with lognormal with a dispersion 
of ~ 0.3. A lognormal distribution of column densi- 
ties naturally ar ises from simulations of turbulence in an 
isothermal ISM dVazquez-Semadenilll994HScalo et al.lll998l; 
iNordlund & Padoanlll999HOstriker et alj|2001l) . with a width 
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FIG. 1. — Star formation and metal enrichment histories for the two syn- 
thetic galaxies considered herein: a passive and a star-forming (SF) galaxy. 



(6) 



that scales with the Mach number. iFischera et alj (12003b 
attempt to f it the observed attenuation law derived by 
ICalzetti et al.l (120001) with a model that includes a clumpy 
ISM. Their best-fit values for the dispersion range from 0.6 < 
a < 2.2. 

Motivated by these results, we adopt a lognormal PDF for 
the distribution of column densities (i.e. P2 in Equation is 
a Gaussian in the variable log(n)). Equation[5]implies a con- 
straint on the relation between the mean and variance of a 
lognormal PDF that we fit with a fifth-order polynomial: 

H w 0.0709 CTdump - 1 .68 of lump + 1 .56 cr;? lump - 

1 ■ 96 Odump + ■ 8 8 6 O'clump > 

where fi = /ii og „ is the logarithmic mean and err lump is the log- 
arithmic variance. This fitting function is valid over the range 
0.0 < fTciump < 1-3. Note that our default assumption of a 
uniform distribution of dust is recovered in the limit where 
Cciump = (i- e - m this limit Pj = 8(n— 1)). The dispersion char- 
acterizing the lognormal PDF is therefore the single variable 
controlling the dumpiness of the ISM in our model. 

In reality one might expect the distribution of column densi- 
ties to be more complex than lognormal. In t he limit that dust 
lies a t the mid-plane of a galactic disk (e.g. Dal canton et alj 
2004), fully one half of the disk stars will experience no at- 
tenuation. For galaxies that contain both a stellar bulge and 
disk, the majority of the bulge stars will also experience lit- 
tle or no attenuation, again under the assumption that most of 
the dust lies at the mid-plane of the disk. While we do not 
include these more realistic g eom etries in our general anal- 
ysis, we will demonstrate in ij4.ll that they can have a very 
important effect on the derived colors of synthetic galaxies. 



2.3. Parameter ranges and priors 

This section describes and motivates the adopted priors on 
each parameter. In all cases we will assume, for simplicity, 
that the prior distribution is fiat between the minimum and 
maximum of the range. 

• Al: Shift in log(Lboi) with respect to the default evo- 
lutionary tracks of TP-AGB stars. There is a lack 
of data on the TP-AGB phase outside of the Galaxy 
and the Magellanic Clouds, making it very difficult 
to even specify a reasonable prior range for this pa- 
rameter. We adopt the following optimisti c estimate: 
-0.2 < Ai < 0.2, and refer the reader to lConrov et alj 
(2009) for motivation of this range. 

• At'. Shift in log(r e ff) with respect to the default evo- 
lutionary tracks of TP-AGB stars. As above, obser- 
vational constraints are lacking for this parameter. We 
adopt the opti mistic range : — 0.1 < A 7- < 0. 1, and refer 
the reader to lConrov et alj (120091) for details. 

• /bhb '■ Fraction of HB stars that are blueward of the red 
clump. It is well known that metal-poor globular clus- 
ters ([Fe/H15 , —1.4) have extended HB morphologies 
dHarrisll 19961) . In addition, there are several less metal- 
poor glob ulars ([Fe/H~|fti - 0.6) with blue horizontal 
branches (Rich et all 119971) . and the most metal-rich 
Gala ctic star cluster kn own, NGC 6791, has /bhb ~ 
0.3 (iKalirai et al.l 120071) . In our Galaxy, several per- 
cent of stars are meta l-poor dZoccali et alj |2003, 2008; 
ISchoerck et al.l 120081) . while in e lliptical galaxies the 
fraction may range from 1 - 10% ( Worthev et al. 199 
iMaraston & Thomas! 120001) . iDorman et alTdl995l 
demonstrated that the observed UV upturn in elliptical 
galaxies can be explained by 5 - 20% blue (extreme) 
horizontal branch fractions. Based on these results, we 
adopt 0.0 < /bhb < 0.2. 

• Sbs- Specific frequency of blue straggler stars, de- 
fined as the number of blue stragglers per unit hori- 
zontal branch star. Typical values for Srs ran ge from 
0. 1 - 1 .0 for globular clusters lPiotto et all ( 120041) . While 
much more challenging to measure, the frequency of 
blue stragglers in the fi eld may be as high as five 
dPreston & Sned en 2000). We adopt the following 
range for Sbs, which will be high if blue stragglers only 
exist in globular clusters, but appropriate if they origi- 
nate from binary star systems, which are ubiquitous in 
the field: < S B s < 2. 

• aiMF- Logarithmic slope of the IMF above IMq. This 
parameter is not included in our total uncertainty bud- 
get, but is inclu ded in Figure 5 to provide intuition for 
it's importance. Rroupa (2001) quote amF = 2.3 ± 0.7 
forM > IMq. The quoted uncertainties are 3a, statis- 
tical. However, as pointed out by Rroupa, the correc- 
tions to a due to unseen binary companions depends 
on a and is of the order of 0.3-0.4. Moreover, direct 
constraints on the IMF come only from the Milky Way 
and Magellanic Clouds, while we require the IMF for 
star clusters of all ages and metallicities. In addition, 
there is increasing indirect evidence that the IMF was 
substantially different when the Universe was younger 
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and in diff erent e nvir onments at low red shift (for a sum- 
mary, see ^2.1.21 and Elmegreen 2009). For these rea- 
sons we adopt the following range: 1.6 < q;imf < 3.0. 

• 6: Power-law index of attenuation cu rve. In a com- 
pariso n to properties of local starbursts, ICharlot & Fall 
(2000) find that the observations can be well-fit with 
0.4 < 5 < 1.3. The upper limit also corresponds to the 
Milky Way extinc tion curve in the optical and near-IR 
( Fitzpatrickl 1 1 9991) . while the lower limit roughly ap- 
proximates the optical region of an extinction curve 
with R = Ay /E(B-V) « 6. Extinction curves with such 
high values of R are seen along d ense sightlines in the 
Galaxy fe.g.. lCardeiri et af]ll98 9\ The attenuation law 
for local starbursts can be characterized b y a power-law 
index of w 0.9 from the UV to w l/im dCalzetti et alj 
2000). Motivated by these considerations, we adopt the 
following range: 0.4 < S < 1.3. 

• f esc : For star clusters, the transition timescale from 
birth clouds to the diffuse ISM. This transition occurs 
either because the birth cloud is destroyed (i.e. pho- 
toevaporated by O-type stars), or because the stars 
wander out of the cloud. The relative importance of 
these two processes is uncertain but will depend on 
the mass of the cloud, in addition to the local star 
formation rate. A variety of observational estimates 
both in the Galaxy and the Local Group suggest that 
this p arameter lies in the range 6.5 < logfe sc /yrs) < 
7.5 dIsraelll978tlWilliarns & McKeell99 7; Blitz & Shul 
ll980tlBritz et alj|2007t iMcKee & Ostrikerll2007l) . This 
range will be adopted herein. 

• n: Optical depth aroun d young stars (where t < 
t esc ). ICharlot & FaTil (120001) found that a variety of data 
from starbu r st gala xies are bracketed by < t\ < 2. 
Humphreys (1978) tabulated V-band optical depths 
for hundreds of supergiants and O-type stars in the 
Galaxy and found a rough range of 0.4 < t\ < 3. 
Observati ons of extragalactic HII regions find simi- 
lar re sults (Isra el & Kennicuttll980l : lvan der Hulst et alj 
1988). We therefore adopt a range of 0.3 < n < 1.5 as 
the range for this parameter. 

• T2'. Optical depth around old stars (where t > f esc ). 
In the context of our dust model, this parameter may 
be considered as a proxy for the dust mass in a galaxy. 
In our total uncertainty budget, this pa rame ter is not al- 
lowed to vary because, as discussed in m.2\ it may, with 
some degree of faith, be predicted in galaxy formation 
models from the surface-averaged column density of 
metals. We vary this parameter in Figure|2]only to pro- 
vide intuition for its importance. 

• c c i ump : Dispersion of the lognormal PDF characteriz- 
ing the distribution of column densities of the ISM. Ow- 
ing to the absence of constraining data beyond the so- 
lar neighborhood, we allow the <7 c i ump to vary between 
0.0 < (T c i ump < 1 .3. The upper limit is set by the practi- 
cal consideration that the effect of cr c i ump on colors tends 
to saturate near this value, and so allowing larger val- 
ues has no effect on observables (see Sj4]for details). A 
value of 0.0 corresponds to a uniform screen of dust. 



3. SYNTHETIC GALAXIES 

As described in the Introduction, our goal is to investigate 
to what extent uncertainties in the SPS technique hamper the 
translation of synthetic galaxies into observables. The out- 
puts of galaxy evolution models typically include star forma- 
tion and metallicity-enrichment histories of galaxies, and, in- 
creasingly, information about the spatial distribution of gas as 
well. 

We make use of the Millennium Simulation database 4 to 
query the outputs from the semi-an alytic galaxy formation 
model of iDe Lu cia & Bl aizotl (120071) . This model couples 
analytic recipes for gas cooling, star formation, supernova 
feedback, AGN heating, galaxy merging and metal enrich- 
ment, to a large, cosmological, dissipa tionless jV-body sim - 
ulation (the 'Millennium Simulation'; Springel et al. 2005). 
The details of their model do not concern us here — we are 
merely interested in extracting a reasonable star formation 
and metal-enrichment history for a typical passive and star- 
forming galaxy. 

Note that spatial information on the scale of a galaxy is not 
available for the galaxy evolution model we consider, in par- 
ticular the spatial distribution of the gas with respect to the 
stars is not a prediction of the model. This is in contrast to 
modern cosmological hydrodynamic simulations. Of course, 
without this detailed spatial information it is extremely diffi- 
cult to accurately account for the attenuation of starlight by 
dust. 

For the bulk of our analysis we will focus on two galaxies 
that represent both passive (non-star forming) and star form- 
ing galaxies. The star formation and metal-enrichment histo- 
ries for these two galaxies are shown in Figure Q] Notice that 
the passive galaxy experienced the bulk of its star formation in 
the distant past (and has experienced no star formation in the 
past w 4 Gyr), while the star-forming galaxy has experienced 
multiple episodes of intense star formation super-imposed on 
a relatively qu iescent level of star formation. Assuming a 
iKroupal d2001) IMF 5 , the passive (star-forming) galaxy has 
a stellar mass at z = of log(M/M ) = 1 1.1(10.9). 

When considering population statistics (i.e luminosity func- 
tions; LFs) we make use of a subsample of the full galaxy 
formation model. We draw galaxies at random from the full 
model with a sampling rate of 10 -4 , resulting in « 7000 syn- 
thetic galaxies. This particular number was chosen as a com- 
promise between adequately sampling the underlying distri- 
bution and computational cost. Since the models do not pre- 
dict dust content, we must specify a recipe for adding dust 
to the galaxy formation model. We make the simple assump- 
tion that T2 = VSFR/2, wher e SFR is measured in M^/yr, and 
then a ssume that ti = r 2 /0.3 (ICharlot & Falll2000l:lKong et alj 
2004). These simple assumpti ons do not deviate dramati - 
cally from the recipes assumed in lDe Lucia & Blaizotl d2007l) . 
where ti is a function of the cold gas mass and metallicity. 
We will be most interested in differential effects, so details of 
the dependence of t 2 on galaxy properties is less of a concern 
here than elsewhere. 



4 http://www.mpa-garching.mpg.de/millennium/ 

5 Given a star formation history, the IMF effects the total stellar mass only 
insofar as it controls the fraction of initial stellar mass formed that is lost due 
to supernova and stellar winds, via the relative proportions of low and high 
mass stars. 



COMPARING GALAXY EVOLUTION MODELS TO OBSERVATIONS 



7 




FIG. 2. — Effec t of d ust parameters on UV, optical, and near-IR colors of a star-forming galaxy (solid lines) and a passive galaxy (dashed lines). The parameters 
are described in £|2.3I The colors are shown as differences from a default model where <5 = 0.7, t\ = 1 .0, Ti = 0.3, ( eS c = 7.0, and cr^u™ = 0.0. In each panel, only 
one parameter is varied. For the passive galaxy, the colors are independent of the parameters T| and f csc because there are no young stars in this galaxy and thus 
these parameters, which only effect dust around young stars, are not relevant. The dependence of the colors on a" c i ump is the same for both galaxies. 



4. RESULTS: UNCERTAINTIES IN TRANSLATING SYNTHETIC 
GALAXIES INTO OBSERVABLES 

This section demonstrates the difficulty of translating the 
synthetic galaxies produced from galaxy evolution models 
into observables. Discussion is focused around three colors, 
one of each sensitive to the ultraviolet, optical, and near- 
infrared regions of the spectrum. We will first investigate 
the dependence of these colors on the parameters in our dust 
model, and then the parameters controlling the SSPs. We will 
then combine all uncertainties to demonstrate their net effect 
on the derived colors, and will explore the effect of these un- 
certainties on model LFs. Discussion is focused on broad - 
band filters available fr om the GALE X dMartin etal.1 120051) . 
SDSS dYork et al.l2000ft . and 2MASS dJarrett et al.ll2000h sur- 
veys. For reference, the filters considered and their effective 
wavelengths are NUV = 2300A, g = 4700A, r = 6200A, and 
K = 2.2/im. 

4.1. Dependence of synthetic galaxy colors on uncertain SPS 

ingredients 

4.1.1. Dust uncertainties 

Figure [2] shows the colors of two synthetic galaxies as a 
function of the parameters of our dust model. For each panel, 
only one parameter is varied. The fiducial set of parameters 



is: 8 = 0.7, n = 1.0, t 2 = 0.3, f esc = 7.0, <7 c i ump = 0.0, A L = 
0.0, A r = 0.0, /bhb = 0.0, S B s = 0.0, and ai MF = 2.3. Several 
important trends are apparent. First, the effect of T2 on the 
colors is the strongest of the parameters explored. This is 
not surprising as t 2 controls the total amount of attenuation 
in the diffuse component. The slope of the attenuation curve, 
5, and the timescale for escape from the birth cloud, f esc , are 
similar in that they have a substantial effect in the UV, with 
increasingly diminishing importance toward the near-IR. The 
parameters t\ and f esc , which control the attenuation around 
young stars, have no effect on the passive galaxy because that 
galaxy has no young stars. The parameter t\ saturates at high 
values for the star-forming galaxy because, as one completely 
extinguishes the light from young systems (where t < f esc ), the 
colors will asymptote to the colors of the older population, 
which in this case are still relatively blue. 

Finally, it is clear that the dumpiness of the ISM, controlled 
by Cciump, has the effect of making colors bluer, and its ef- 
fect tends to saturate at cr c \ ump ~ 1.3. Colors become bluer 
as crciump increases because the average optical depth must de- 
crease for a clumpier medium (see Equation|6]). The effect sat- 
urates when the average optical depth approaches zero, which 
occurs at a c i ump ~ 1.3 (again see Equation|6]l. In this case most 
stars are unobscured and therefore a further increase in cr c i ump 
cannot alter the integrated colors of the whole population. 

The importance of the dumpiness of the ISM is explored 
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FIG. 3. — Dependence of colors on the dispersion in column densities of 
the ISM (where o- c i um p = 0.0 corresponds to a uniform screen of dust). Each 



relation is normalized to the colors at a, 



clump 



= 0.0. The relations are shown 



for several values of the optical depth parameter r 2 , as labeled in the figure. 
The t 2 = 0.3 case is identical to the far-right column of Figure [2] These 
results are shown for the star-forming galaxy only. This figure demonstrates 
that the importance of carefully handling the density fluctuations in the ISM 
increases dramatically at higher mean optical depths, i.e. for larger total dust 



further in Figure [3] This figure is similar to the far-right col- 
umn of Figure [21 except here we show results only for the 
star-forming galaxy, for varying values of T2 (notice also that 
the range of the y-axes are different). It is clear that modest 
increases in result in a substantially stronger dependence of 
colors on the dumpiness of the ISM. This trend is due to the 
fact that increasing u c i ump decreases the effective dust opacity 
to approximately it's limiting value of zero (as can be seen by 
comparing columns three and five in Figure |2). Since varia- 
tion of Uciump can be interpreted as a variation of the effective 
dust opacity from ti to w 0.0, increasing T2 will clearly in- 
crease the effect of a c i ump on the derived colors. 

The implications of Figure [3] should not be underesti- 
mated. For galaxies that are moderately dusty {tx > 1), 
the geometry of the dust is enormously important when 
translating physical properties into observations, and visa 
versa. Note also that the dependence of colors on cr c i ump 
is strongest for modest deviations from a uniform distribu- 
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FIG. 4. — Dependence of colors on the fraction of starlight unobscured by 
dust, / (a fraction of 0.0 corresponds to our standard model where all starlight 
is subject to attenuation). Each relation is normalized to the colors at / = 0.0. 
The relations are shown for three values of the optical depth parameter r%, 
as labeled in the figure. These results are shown for the star-forming galaxy 
only. 

tion of dust. Observations suggest that in the Milky Way 
Cciump ~ 0.3 ((Berkhuiisen & Fletcher 2008), where the de- 
pendence of color on <7 c i U mp is strongest. This suggests that 
the detailed distribution of dust with respect to the stars must 
be known to high precision. Given this sensitivity, it is rea- 
sonable to suspect that deviations from our assumption of a 
lognormal distribution of column densities will also be im- 
portant. 

In Figure|4]we explore further the importance of the large- 
scale spatial distribution of the dust with respect to the stars. 
In this figure we modify our default dust model by allowing 
for a fraction of starlight to be unobscured by diffuse dust (the 
aspect of the dust model regarding young stars is unchanged). 
As discussed in £12.2.11 realistic galaxies will likely have a 
substantial fraction of starlight unobscured by dust. Consider 
again the example of a galaxy composed of a disk and a bulge 
with a majority of its cold gas and dust at the midplane of 
the disk. At any viewing angle, roughly one half of the disk 
stars, and the majority of the bulge stars, will be unobscured 
by dust. Unobscured light fractions as high as 50% will there- 
fore be common in galaxies containing dust. Figure [4] fur- 
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FIG. 5. — Effec t of S SP parameters on UV, optical, and near-IR colors of a star-forming galaxy (solid lines) and a passive galaxy (dashed lines). The parameters 
are described in £|2.3I The colors are shown as differences from a default model with Al = 0.0, At = 0.0, /bhb = 0.0, Sbs = 0.0, oimf = 2.3, and with the default 
dust parameters as in Figure[2] The population of blue horizontal branch stars and blue straggler stars, controlled by the parameters /bhb and Sbs, have little 
effect on the colors of the star-forming galaxy because this galaxy contains a large number of young hot stars that outshine these hot, evolved populations. The 
IMF slope, aiMF, is only varied for masses > IMq. 



ther demonstrates that fruitfully comparing galaxy models to 
observations requires a detailed understanding of the relation 
between dust and stars. 



4.1.2. Stellar evolution and IMF uncertainties 

Figure [5] shows the colors of two synthetic galaxies as a 
function of the parameters controlling the SSPs, i.e. param- 
eters governing late stages of stellar evolution and the loga- 
rithmic slope of the IMF at high masses (> IMq). It is clear 
that the TP-AGB parameters, Al and At, have little or no 
impact in the UV and optical, but a large effect in the near- 
IR. The effect is larger for the star-forming galaxy because 
this galaxy has a larger fraction of intermediate age stars, and 
the TP-AGB phase contributes substantially to the bolometric 
luminosity of inter mediate age populations (Maraston 2005; 
IConroy et alJl2009h . 

The morphology of the horizontal branch and the frequency 
of blue straggler stars, parameterized by /bhb and Sbs, respec- 
tively, only impact the colors of the passive galaxy. These pa- 
rameters have no effect on the star-forming galaxy because 
this galaxy has young, hot stars that outshine the hot, but 
less luminous, blue straggler and BHB stars. The impact 
of these parameters on the UV light for the passive galaxy 
is quite strong. Even a modest number of blue horizontal 
branch or blue straggler stars can produce changes in the UV 
of several tenths of a magnitude. The impact of these exotic 



stars on the UV spectrum of passive galaxies has been dis- 
cussed for decades as a possible explan ation for the UV up- 
turn phenomenon in elliptical gala xies (Bur stein et al.ll 1988k 
Doi- man et alJ[l995l:lHan et alj|200l . 

The final column of Figure [5] shows the effect of the log- 
arithmic slope of the IMF on the derived colors of synthetic 
galaxies. The slope is only allowed to vary for stellar masses 
> 1M Q . The explanation of the trends seen in the figure are 
qualitatively different for the two galaxy types. The trends 
seen for the passive galaxy are exp lained by inspectio n of the 
IMF-dependence of the SSPs (see lConrov et alj|2009l for de- 
tails). For wavelengths < 5000A, the SED is sensitive to stars 
at a particular mass (i.e. stars near the main-sequence turn- 
off mass). However, at wavelengths > 5000A, the SED is sen- 
sitive to both the red giant branch (RGB) and asymptotic giant 
branches (AGBs), and in particular the TP-AGB phase. For 
a coeval set of stars, these two phases of stellar evolution are 
inhabited by stars of slightly different masses (stars along the 
AGB are more massive than stars along the RGB), and there- 
fore the relative weights given to these phases will depend on 
the IMF. In essence, a steeper IMF favors the RGB over the 
AGB, and since the AGB is at lower effective temperature, a 
steeper IMF results in slightly bluer near-IR colors. 

The trends of «imf with color for the star-forming galaxy 
are both much stronger and, in the optical and near-IR, in the 
opposite sense compared to the passive galaxy. These trends 
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FIG. 6. — Distribution of A^-band magnitudes and broad-band colors arising from marginalizing over uncertainties in the dust model and stellar evolution. All 
quantities are plotted as differences from our default model. These distributions are created by marginalizing over all parameters shown in Figures|2]and|5] except 
for ctimf and T2, for reasons discussed in the text. Left Panel: Distributions for a star-forming galaxy. Results are shown for two values of ti: T2 = 0.3 (solid 
lines), and T2 = 1.0 (dashed lines). The difference between the two lines is driven by the dumpiness of the ISM, and is explained by the results in Figure|3] Right 
Panel: Uncertainties for a passive galaxy. Results are shown for two values of r^. T2 = 0.3 (solid lines), and t-i = 0.0 (dashed lines). Note that in the latter case 
this implies that the galaxy is entirely dust-free. 



are due to the fact that the colors of a star-forming galaxy are 
sensitive to stars of a range of masses. This can be understood 
from the following example. Imagine that a galaxy consists 
of two SSPs: an old and a young population. As mentioned 
above, the light from each component is dominated by stars 
of a given (different) mass, and so the relative importance of 
each young and old component is determined by both the star 
formation history, which specifies the fraction of mass formed 
in each component, and the IMF, which specifies the relative 
weights given to the mass intervals that dominate each com- 
ponent. One can see this another way by considering an ana- 
lytic representation of the combined flux: 

F\ ~ $(M 1 )*(/ 1 )S A (M 1 ,? 1 ) + $(M 2 )*fe)5 A (M 2) ? 2 ), (7) 

which is qualitatively analogous to Equation [Tj without dust. 
This equation is approximate insofar as we are assuming that 
the old population, born at time t\, is dominated by stars of 
mass Mi, while the young population born at time f 2 is domi- 
nated by stars of mass M 2 . The logarithmic slope of the IMF 
determines the ratio $(M2)/$(Mi), and it is therefore clear 
that the IMF directly effects the emergent flux for a fixed star 
formation history. 

In the context of the above example, increasing ums results 
in a redder SED because the young (blue) component is sensi- 
tive to high-mass stars, while the old (red) component is sen- 
sitive to low-mass stars, and a steeper IMF favors low-mass 
over high-mass stars. This simple example serves to qualita- 
tively explain the trends observed in the far-right column of 
Figure [5] 



4.2. Total uncertainties in UV, optical, and near-IR colors 

In the previous section we explored the dependence of syn- 
thetic galaxy colors on various uncertain aspects of the SPS 
technique, including uncertainties in the dust model, stellar 
evolution, and the IMF. We are now in a position to investigate 
the combined effects of these uncertainties on the broad-band 
colors of synthetic galaxies. 

Figure|6]shows the distribution of colors and /T-band mag- 
nitudes for two synthetic galaxies. The distributions were ob- 
tained by marginalizing over the uncertain aspects of SPS dis- 
cussed in the previous section. All quantities are plotted as 
differences with respect to our default model. We marginal- 
ize over all parameters shown in Figures [2] and [5] except for 
aiMF and r 2 . We choose to fix r 2 because, if we interpret this 
parameter as probing the total dust mass, then we may hope 
that galaxy evolution models are able to predict this quantity, 
if not at present then in the future. In other words, our aim in 
this work is to marginalize over aspects of SPS that must be 
assumed when translating models into observations. If one or 
more parameters discussed in the previous section is readily 
predicted by such models, then such parameters should obvi- 
ously not be marginalized over. 

We choose to fix oimf simply because its effect on the re- 
sulting colors is so large that it would dominate the error bud- 
get in all cases. The uncertainties in broad-band colors in- 
ferred from Figure [6] can thus be safely interpreted as lower 
limits to the true uncertainty. 

Figure [6] shows results for different values of r 2 . For both 
galaxies we consider our default value of r 2 = 0.3. For the 
star-forming galaxy we also consider r 2 = 1 .0, while for the 
passive galaxy we consider r 2 = 0.0. We vary r 2 for the star- 
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forming galaxy to demonstrate the fact that the relationship 
between cr c iumn and the resulting colors itself depends on T2 
(cf. Figure |3). Larger values of r 2 imply that variation of 
Cciump will have a progressively larger impact on the derived 
colors, as is evident in the broader distributions in Figure [6] 
forr 2 = 1.0. 

Passive galaxies are often observed to contain little dust 
(t 2 < 0.1; see e.g iGoudfrooii et alj|1994t iFerrari et al.lll999t 
iTemi etalJ 2004; D raine et all 120071) . We thus also include 
T2 = 0.0 in Figure|6]for the passive galaxy, which implies that 
the parameters cr c i urnp and 5 have no effect in the marginal- 
ization. The removal of dust decreases the uncertainties only 
slightly. This is because, for the passive galaxy, uncertainties 
are dominated by the uncertainties in stellar evolution. 

LFs constructed from a subsample of the full galaxy for- 
mation model are shown in Figure [7] both in the /f-band 
and M/V-band. The LFs are computed from the model star 
formation and metallicity histories, using the SPS model de- 
scribed in $2] LFs are constructed for a variety of assump- 
tions regarding the stellar evolution and dust model parame- 
ters. Note that the LFs were constructed from the exact same 
galaxy sample in each case, so the variations seen in Figure 
are entirely due to the different modeling assumptions (as 
opposed to Poisson or sample variance). 

Clearly, the TP-AGB parameters Ai and A7- have a large 
and non-trivial effect on the Zf-band LF. The dust parameters 
have a negligible effect on the /f-band LF, and so their vari- 
ation is not included in the /f-band LF. In contrast, the dust 
model parameters have a dramatic effect on the NUV-band 
LF. The TP-AGB parameters have no effect in the NUV. 
While not shown, sensitivity to the ensemble of parameters 
is roughly minimized in the r-band, suggesting that the opti- 
cal portion of the SED is the most robust to SPS uncertainties. 

Notice that for both LFs in Figure [7] the shape of the LF 
changes for different model parameters. This is due to the 
fact that the effects of these parameters vary with stellar pop- 
ulation age. The effects of SPS uncertainties therefore induce 
complex systematic uncertainties in model LFs in the UV and 
near-IR that must be carefully handled if robust constraints on 
galaxy formation models are to be achieved. 

5. DISCUSSION 
5.1. How do we interpret these uncertainties? 

We have demonstrated that the broad-band colors of syn- 
thetic galaxies carry significant uncertainties due to uncertain 
aspects of SPS. These uncertainties must thus be incorporated 
into any analysis that attempts to compare models to observa- 
tions. But how are these uncertainties to be interpreted? 

Unfortunately, it is not clear whether or not one should treat 
these uncertainties as statistical or systematic. By statisti- 
cal we mean that the uncertain SPS parameter values vary 
stochastically from galaxy to galaxy, whereas by systematic 
we mean that these parameters vary systematically across 
galaxy types. Uncertainties associated with the dust model, 
including the attenuation curve and large-scale geometry of 
the ISM, can probably be interpreted as statistical uncertain- 
ties. These aspects depend in very complex and uncertain 
ways on the geometry of the stars with respect to the dust, 
the local population of O-type stars, the gas-phase metallic- 
ity, and other physical properties of the ISM. 

Uncertainties associated with stellar evolution are more 
likely systematic, although there may be a statistical compo- 
nent. One might imagine that the fraction of blue straggler 



or horizontal branch stars, or the evolution of TP-AGB stars, 
would not vary stochastically from galaxy to galaxy since 
stellar evolution is believed to depend on a small number of 
variables. However, even within globular clusters there ap- 
pears to be a need for stochastic mass-loss along the RGB 
in order to explain the observed broa d range in tempera- 
tures along the ho rizontal branch (e.g.. iDorman et allll995l ; 
iKalirai et al.ll2007h . Thus, in addition to the systematic varia- 
tion of these post-main sequence stars with global quantities 
such as metallicity and age, one might indeed expect a statis- 
tical uncertainty as well. 

Uncertainties in the IMF may be plausibly be either statis- 
tical or systematic, depending on the dominant mechanisms 
responsible for the shape of the IMF (i.e., is gas-phase metal- 
licity or ISM pressure more important?). 

Whether or not the resulting color uncertainties are predom- 
inantly statistical or systematic obviously has very different 
implications for their interpretation. For the statistical errors, 
the interpretation is straightforward. Systematic uncertainties 
are rather more troubling in this context because systematics 
will be correlated with galaxy properties such as age, metal- 
licity, and/or mass. For example, if future observations of 
TP-AGB stars imply that A^ = 0.3, then the colors will shift 
in an age-dependent (and because mass correlates with age, 
a mass-dependent) way. Thus, systematic uncertainties will 
tend to introduce not only overall shifts but also tilts in ob- 
servational quantities derived from models such as the color- 
magnitude diagram and LFs. 

With the preceding discussion in mind, we return to Fig- 
ure [6] It should be clear that the distributions shown in this 
figure are not to be interpreted as PDFs for colors given a 
star formation history. Consider a simple example. If we 
had chosen to marginalize over nine parameters that had al- 
most no effect on a particular color and one parameter that 
had a large effect, then the distribution in colors would be 
very strongly peaked at A(color) w 0, with broad wings at 
low frequency, because each parameter is sampled equally. 
Even with a strongly peaked distribution, one would not con- 
clude that the uncertainty on the color is small, since that one 
important parameter produces such large uncertainties. The 
implication is that the only robust way to interpret these dis- 
tributions is to consider their entire width as a simple measure 
of the resolution with which one may accurately predict colors 
given star formation and metal enrichment histories. 

5.2. Why the IMF matters 

It is often assumed that results are insensitive to the IMF 
so long as both modelers and observers use the same IMF to 
translate between fluxes and physical parameters. This as- 
sumption is incorrect. Galaxy evolution models predict star 
formation histories (SFHs) that are independent of the IMF. 
Upon integration, these SFHs provide total stellar masses for 
model galaxies. In this context, the IMF enters only in the es- 
timation of the fraction of mass returned to the ISM via winds 
and supernovae. 

In contrast, since the IMF has a direct effect on the colors 
of star-forming galaxies, it will effect the SFRs and stellar 
masses inferred from observations, and will effect modelers' 
ability to translate and compare their models to observations. 
It is important to realize that this is true even if the low-mass 
end of the IMF — which has a negligible effect on the inte- 
grated light from galaxies — is fixed. It is the slope of the IMF 
at masses > 1M Q that most affects the colors of star-forming 
galaxies. 
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FIG. 7. — LFs of synthetic galaxies from the galaxy formation model of De Lucia & Blaizot 1 2007). SFHs and metallicity histories are provided by the model. 
For the same set of galaxies, LFs are constructed in the K— band and NUV— band for a variety of assumptions regarding the dust model (parameterized by 8, 
fclumpi ar| d fesc), importance of the TP-AGB stellar evolutionary phase (A/, and At), and fraction of blue HB stars (/bhb)- It is clear that for the same galaxy 
formation model, the predicted LFs can span a wide range depending on the details of how one translates the model galaxies into the observational plane. 



The IMF is thus extraordinarily important when attempting 
to compare models to observations. As we have demonstrated 
herein, the measured uncertainty in the IMF in the solar neigh- 
borhood is sufficiently large to adversely effect comparisons 
between models and observations. A more pessimistic assess- 
ment of our knowledge of the IMF in external galaxies (and 
at earlier epochs), will only exacerbate this problem. 

5.3. What is to be done? 

Despite the magnitude of the uncertainties presented herein, 
there is, we believe, a clear way to proceed in the comparison 
between models and observations. 

As emphasized in the present work, the principle difficulty 
in translating models into observables is the array of important 
aspects of SPS that are neither well-constrained observation- 
ally nor well-understood theoretically. We are thus forced to 
adopt rough ranges for these uncertainties — rough because 
uncertainties are difficult to quantify. However, this does not 
imply that each observed galaxy can be equally well-fit by 
the full range of parameters. In other words, it is precisely 
because these parameters impact the broad-band colors of 
galaxies that they may be constrained, to some extent, directly 
by the data. 

We therefore advocate using the SPS technique to estimate 
the basic physical properties of observed galaxies, including 
star formation rates and stellar masses. The SPS technique 
must include a marginalization over the uncertain aspects dis- 
cussed herein. The derived properties will therefore be robust, 
though they will carry larger errors. Such an approach re- 
quires a pan-chromatic view of a large number of galaxies at 
multiple epochs. The physical properties of observed galaxies 
can then be directly compared to galaxy evolution models, so 
long as proper attention is given to the impact of the associ- 
ated uncertainties. We believe that this approach will afford 



a more accurate and stringent constraint on models than the 
approach of converting models into observables. 

This recommendation is guided largely by the manner in 
which modelers, in practice, compare their predictions to ob- 
servations. Models are often tested against observed LFs, 
color-magnitude diagrams, and luminosity-dependent clus- 
tering. We instead suggest that models should be tested 
against observationally-constrained physical relations includ- 
ing mass functions, mass-dependent clustering, and SFR- 
mass correlations. As these observational results have become 
available in the past few years, modelers have begun to use 
them as constraints. We believe that this is the better way to 
proceed, as long as two conditions are met: 1) the observa- 
tional results contain an accurate accounting of all the rele- 
vant uncertainties, including the correlations between derived 
products, and 2) the comparison between model and obser- 
vations takes careful account of the resulting uncertainties on 
the derived physical properties. 

Proponents of the alternative approach, namely that mod- 
elers translate their results into the observational plane, may 
counter that in the models the SFH and metallicity of the stars 
is known precisely, whereas in the data these aspects must 
be estimated fr om the data. The age-metallicity degeneracy 
(IWortheyll9 94). for example, may present an ambiguity in the 
fitted data, whereas none would arise by converting the mod- 
els into the observational plane. We believe that the fitting the 
SFH and metallicity in the data, in addition to fitting for dust 
and SPS uncertainties, is still the preferred approach, rather 
than marginalizing over the uncertain dust model parameters 
and SPS uncertainties when attempting to robustly translate 
model results into observation al predictions. This v iew is sup- 
ported by the recent work of Galla zzT& Belli (120091) who find 
that in the absence of dust, uncertainties in fitting SFHs and 
metallicities translate into a relatively minor uncertainty in the 
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stellar mass-to-light ratio. These authors conclude that dust 
obscuration and SPS uncertainties remain the dominant un- 
certainties when estimating physical parameters of galaxies. 

6. SUMMARY 

We have investigated some of the uncertainties associated 
with translating synthetic galaxies into observables, focus- 
ing on broad-band UV through near-IR photometry available 
from the GALEX, SDSS, and 2MASS surveys. The uncer- 
tainties associated with stellar evolution and dust are impor- 
tant and effect passive and star-forming galaxies in different 
ways. 

For star-forming galaxies, both uncertainties in dust and in 
the TP-AGB phase of stellar evolution result in substantial 
uncertainties in UV, optical, and near-IR colors. Uncertain- 
ties in the spatial distribution of dust with respect to stars play 
an especially important role, as do uncertainties in the dust 
attenuation law. 

The uncertainties in the colors of passive galaxies are dom- 
inated by uncertainties in stellar evolution, including the TP- 
AGB phase, blue stragglers, and the morphology of the hor- 
izontal branch. The uncertainties in the dust model have a 
minor effect because the total dust content of passi ve galaxies 
is observed to be low, though non-zero (e.g .[Goudfrooii et al. 
[l^^lFerrari"etal] [r999: Drai ne et alj|200% . 

We have also explored the impact of the IMF on the derived 
colors of galaxies, although we have not included this source 
of uncertainty in the total error budget. Varying the logarith- 



mic slope of the IMF has different effects for passive and star- 
forming galaxies. The magnitude of the effect is much larger 
for star-forming galaxies, and needs to be taken seriously by 
both modelers and observers. The importance of the IMF in 
this re gard has been known since at least the work of lTinslevI 
(1980) — we have explored it here primarily to (re)call at- 
tention to its importance when comparing models to observa- 
tions. 

The ranges of uncertainties investigated here are largely op- 
timistic (i.e. they are relatively small), and yet their influence 
on our ability to fruitfully compare models to observations 
is substantial. Many of these uncertainties are systematic, 
further complicating their interpretation. The comparison of 
models to observations can be done with confidence only once 
these uncertainties are carefully included in the analysis. The 
results from this work highlight the urgent need for obser- 
vations capable of diminishing the various uncertainties dis- 
cussed herein. 
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